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Abstract 

We present the path integral representation of the generating function for classical exclusive 
particle systems. By introducing hard-core bosonic creation and annihilation operators and ap- 
propriate commutation relations, we construct the Fock space structure. Using the state vector, 
the generating function is defined and the master equation of the system is transformed into the 
equation for the generating function. Finally, the solution of the linear equation for the generating 
function is derived in the form of the path integral. Applying the formalism, the equivalence of 
reaction-diffusion processes of single species and two species is described. 
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I. INTRODUCTION 



Recently, nonequilibrium systems such as reaction-diffusion systems, driven lattice sys- 
tems, and cellular automata have been intensely investigated, not only because of their con- 
nections to a variety of important physical phenomena (nonequilibrium phase transitions, 
long-range correlations, and self-organized criticality), but also because of the analytic chal- 
lenge due to the lack of a general formalism to describe nonequilibrium systems [1] . Although 
numerical approaches have played a major role in the investigation of nonequilibrium sys- 
tems and have been successful in many areas, a general field-theoretic formalism has been 
constantly sought because, along with a renormahzat ion-group analysis, it would provide 
a systematic tool to evaluate physical observables in the scaling regime Q]. It is difficult 
to investigate nonequilibrium systems analytically at the microscopic scale because there 
are too many microscopic variables involved, and microscopic variables take discrete values. 
A field-theoretic strategy is to find, starting from the microscopic rules, the equations of 
motion of well-chosen coarse-grained variables. These mesoscopic variables take continuous 
values and vary continuously in space and time. They describe the system at the meso- 
scopic scale, while keeping track of all the fiuctuations related to the microscopic degrees of 
freedom. Doi first introduced a field-theoretic formalism for react ion- diffusion systems using 
the bosonic coherent state path integral P|, and several authors revived the formalism and 
incorporated renormalization-group approaches to the description of the anomalous kinetics 
in reaction-diffusion systems and the stochastic sandpile model [J, , 0, Q, 0, 0] • 

Despite the success achieved by the bosonic field theory for reaction-diffusion systems, 
some nonequilibrium systems with exclusive particles cannot be analyzed correctly by the 
bosonic field theory. Driven reaction-diffusion systems jiol . 11 1, multispecies adsorption 
models [l^l, and driven lattice gases (DLG) [l^ are some examples to which the bosonic 
field theory cannot be applied. In these systems, the exclusive property of the particles is 
important. Several attempts have been made to incorporate the exclusive property into a 
field theory. One way to take the exclusive property into account is to exploit the physical 
knowledge of the system phenomenologically as Zia and Schmittmann have done for the 
DLG. This method is an approximation and cannot be systematic. Other approaches pro- 
posed by Brunei et al. 14j and Bares and Mobilia are formulated using the fermionic 
path integral and are rather difficult to analyze and to extend for higher spatial dimen- 
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sions or multispecies processes. Recently, van Wijland extended the bosonic field theory by 
introducing the exclusion constraint operator to take the exclusive property into account 
Q]. He used bosonic creation and annihilation operators to construct the path integral by 
considering the exclusive property using the exclusion constraint operator. 

The authors previously presented a hard-core bosonic field theory using hard-core bosonic 



creation and annihilation operators [17| . It has been successfully applied to the asymmetric 
exclusion process and several restricted solid-on-solid-type growth models to provide the 
correct Langevin-type equations of motion and the proper path integral formulation jisl ]. 
However, the formalism assumed the existence of the Fokker-Planck equations for processes 
a priori and derived the Langevin type equations of motion. In this paper, we develop a 
more general path integral formalism for nonequilibrium systems with exclusive particles. 
In what follows, we present the derivation of the path integral formulation for systems 
with exclusive particles using the generating function of the correlation functions, and we 
will illustrate how the formalism can be applied to various reaction-diffusion systems to 
establish nonequilibrium universality classes. 

II. FORMALISM 

The dynamics of a stochastic system are usually described by the master equation gov- 
erning the time evolution of the probability P{C;t). P{C;t) is the probability of a system 



being in a microstate C at time t. The master equation takes the form lia |20|, |21| 

d_ 

dt 



^P{C- t) = J2 {^cc'PiC; t) - WacP{C- 1)} , (1) 



CYC 

where Wac is the transition rate, or transition probability per unit time, from microstate C 
to C . Introducing the matrix Hqci 



Hcc = Sec ^ Wa'c - Wcc, (2) 

C" 

Eq. ((T)) can be written succinctly as 

^P{C-t) = -Y,Hcc'P{C'-t). (3) 

c 

Due to the resemblance of the master equation to the Schrodinger equation, it is convenient 
to introduce Dirac's bra-ket notation and an orthonormal basis formed by microstates 



Introducing the state vector \P;t) = P{C; t)\C) and the Hamiltonian operator H whose 
component is Hqc in this orthonormal basis, the master equation can be written in the form 
of an imaginary time Schrodinger equation, 

^^\P;t)^-H\P;t). (4) 

In this representation, the average of a physical quantity can be written as 

{0) = J20{C)P{C;t) = {-\d\P;t), (5) 
c 

where O is an operator with the elements {C'\0\C) — Scfi'0{C) and the projection state (-1 
is defined as 

c 

A natural way to deal with a stochastic system in which the particle number at each site 
varies is to introduce a Fock-space-hke structure. When a system has M different particle 
species with the exclusive property, at each site there are (M+ 1) possible states: a vacuum 
state and M occupied states depending on particle species. Using Greek indices for the 
kinds of particle species — a vacuum state is regarded as a new species and is reserved for 
a vacuum state — and bold Latin indices for site locations, wc denote a state at the site n 
by \an), where an is the index of species occupying the site n and goes from to M. For 
microstates of the whole system, we work in a phase space which is composed of the direct 
product of single-site microstates of all sites {n}, 

\{N})^(^M. (7) 

n 

At each site n, creation and annihilation of the exclusive particles are described by creation 
and annihilation operators as follows: 

aa,n^ \On){an\, al^n ^ \o^n) {On\ («7^0), (8) 

with the closure relation 

M 

X] |aTi)(Qiri| = / (9) 

a=0 

and the orthonormality condition 

{an\Pn)=Sap, Q;,/3 = 0,1,...,M. (10) 



These operators obey the hard-core bosonic commutation relations 

aa,nai3,n = oj, ^ajj^^ = 0, (Ua) 

M 
7=1 

[aa,ri, ajj.rn,] = [«a,ri, a/3,rTi] = {u ^ m) , (llc) 

where Eqs. (|lla|) and ()llb|) prohibit double occupancy at a single site even for different 
species — hence the nomenclature "hard core." Equation pic|) states that any two operators 
defined at different sites commute — hence the name "boson." As usual, Na,n = (ia,n^a,n is 
the number operator whose eigenvalues are (vacant) and 1 (occupied). 

To each state of the system, we can associate the state vector |P; t), considered an element 
of Fock space, 

|P;t) = 5^P({iV};t)|{iV}), (12) 

m 

where |{A^}) is the microstate defined in Eq. (jZj). The microstate can be written in terms 
of hard-core boson operators as 

M 

\{N}) = UWUr-m), 

..=1 (13) 

(wi=({o}inn(^->-)''"-' 

n a=l 

where Na^n is the eigenvalue of the number operator corresponding to the eigenstate {N}) 
{Na^n\{N}) = Na^n\{N})). Usiug hard-core boson operators, the projection state can be 
written as 

(■i=({o}infi+EO- (14) 



a=l 

Now we introduce the generating function F{{ip}]t), defined as 



Fm;t) = {mP;t), (15) 

where 

/ M 
(Ml=(8) (0.|+$^(«r^|/(^a,, 



a=l 
M 



(16) 



a=l 



and / is an analytic function of its argument. F is the generating function in that by 
differentiating F with respect to / and setting / = 1, we can find all kinds of correlation 
functions of the particle density. For example, 

d 



|P;t). 



(17) 



The exact form of / does not matter. In what follows, we mainly use three kinds of functions. 



f{v) = (18a) 

f{v) = 1 + (18b) 

/(^) = exp((^), (18c) 

and these will be used in different contexts. 

By differentiating Eq. (fT3j) with respect to t, we obtain 

^^Fm-t) = -m\H\P-t). (19) 

Since our main goal is to find a path integral representation of -F({(^};t) from Eq. (fT!]|l . 

we need to find a linear partial differential equation about F({(^};t). This can be achieved 
using the following properties of ({<^}|: 

mValr. = f{^c.,r.)Mm. (20a) 

{m<r. = ...^ . (Ml, (20b) 

d 

/(^a,r.) ,,,_ (20c) 

{m<r.a\^r^ = ^a^Mm. (20d) 

where 

(^"^ShT^ 

is the projection operator to the vacuum state at site n. Relations (pUj) yield the partial 
differential equation for the generating function. 



|F(M;t) = -C (m, 1^1) F(M;t), (21) 

where C takes the normal-ordered form, that is, all <^'s are located to the left side of any 
d/dtp. We call C an evolution operator. Since this is a linear equation, we can write the 



path integral solution of F{{(p};t). The path integral solution of Eq. (|7T|) with any of the 
prescription Eqs. (fTH|) is 

dr]a,ndf] 



with 



where 



and 



Fm;t)= /n 



s 



-sm,{'p},t)+m-wm+{vHm} 



dt 



d 



XI fa,r^-^^a,r^ + C{{(p{t)}, {(fit)}) 



(22) 



(23) 



(24) 



(25) 



Taking a continuum limit and keeping the most relevant terms, we arrive at mesoscopic 
action which is equivalent to the microscopic master equation. The long-time properties are 
extracted by studying the action using renormalization-group theory. 



III. EQUIVALENCE BETWEEN STOCHASTIC SYSTEMS 

The equivalence between stochastic systems has been usually studied using similarity 
transformation for single-species react ion- diffusion systems [22]. This section shows the 
equivalence between stochastic systems using the evolution equation of the generating func- 
tion in Eq. (PT|) with the prescription Eq. ()18b|) instead of using the similarity transfor- 
mation. The key mechanism to deduce the equivalence between stochastic systems is the 
rescaling of the field (p (see below), which is simple enough to be applicable to higher- 
dimensional systems and multispecies systems. 

To begin with, we consider a single-species reaction-diffusion model which is defined on 
a d-dimensional hypercubic lattice with diffusion, pair annihilation, coalescence, death, and 
single-particle branching. Particles move with a diffusion constant D. When two particles 
form a nearest-neighbor pair, both of them are annihilated with rate A or one of them is 
removed with rate 1]. Additionally, a single particle is annihilated spontaneously with rate 
6 and a particle- vacant pair becomes a particle-particle pair with rate a. The dynamics 
are summarized in Table lU When 6 = a = 0, this model corresponds to the single-species 
annihilation and coalescence model. The model with D = X = ri = Ois the well-known 
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contact process 23|. When all processes are present, this model is known to show an 
absorbing phase transition which shows the same critical behavior as that of the directed 



24|. 



percolation 

The Hamiltonian of this reaction-diffusion system is H = "^ni^n ~^ ^n~^ ^n~^ ^n~^ ^n) 
with 

d 

Hn= [pri(l - Pri+eJ " a^a^n+e, + (1 - Pn.) Pn.+e, " al^Ori+eJ : 



i=l 
d 

i=l 
d 



i=l 
d 

Hn=7^Yl l.P''^^ ~ ~ Pri^l+e, + (1 - Pn)pn+e, - ^IiPri+eJ : 



(26a) 

(26b) 

(26c) 
(26d) 
(26e) 



1=1 



where pn = CLlfin is the number operator and we drop the unnecessary species index. The 
operators satisfy the commutation relations Eqs. (Illj) with M = 1. Following the procedure 
explained in Sec. IHl we can find the evolution operator C with the prescription Eq. ()18b|) . 
If we write £ as £ = X]f=i ^^n^i-, we obtain 

£ri,i = {'^D + (5)<^„— ( ^ + 77 ) ( VnTr^ \- <^ri+e.77— 



+ {\ + ri + a - 2D)Lpn(^ri+e,— 



[ifn + ^n+e 



92 



(27) 



d d 

+ 



If we use an uncorrelated initial condition with density Pq, the initial-state vector and the 
initial generating function can be written as 



|P;0) = n[(l-Po)+Poay |0), 

n 

^(M;0) = n(i + Po^«). 



(28) 



Hence Eqs. flTTf) and (PHj) along with Eq. (PT|) fully specify the above reaction-diffusion 
system. 
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Let us assume that the solution of Eq. (PTj) is written as F{{(p}, D, A, i], S, a, po; t). Rescal- 
ing the field ip = iiip' (// > 0), F is modified to 

D, X, 7], S, a, Po; t) = F{{^'}, D, A, fj, ~6, a, ppo; t), (29) 

where the relations between parameters with and without a tilde are found by setting ip = 
pip' in Eq. (P7j). which read 

a = per, 
2D + 6 = 2D + S, 

^ + f = ^ + |, (30) 

A+?+?=irA+?+" 



2 2 p\ 2 2 
X + fi + a-2D = X + r] + a-2D. 

Since F on the right-hand side of Eq. can be regarded as a solution of Eq. (PT|) with 
tilded rates and initial density ppo, the two systems connected by Eq. share the same 
generating function. The relation of the correlation function can be found by differentiating 
F with (fi, which reads 

Ck{{x}, D, A, 1], 5, a, Po; t) = p'''Cki{x}, D, A, r/, S, a, ppo; t), (31) 

where Ck (Ck) is the /c-point correlation functions of the reaction-diffusion systems with the 
untilded (tilded) transition rates. Thus, for arbitrary p, which ensures all tilded parameters 
are non-negative, we can find the equivalent stochastic systems to the system with untilded 
parameters. 

Let us find the equivalent systems to the single-species pair annihilation model [rj = a = 
S = 0). For given p, the tilded rates are found as 

D = D,a = 0,5 = 0,fi = 2 ( 1 - - ) A,A = ( - - 1 ) A, (32) 



p J \p 

where 1 < /i < 2 should be satisfied to have the physical meaning. If we choose p = 2 
in Eq. (j32p . the tilded rates become D = D, fj = X, and 5 = 5" = A = 0, which are the 
transition rates of the single-species coalescence model with the initial density po = 2po 
(of course po should not be larger than i). Hence all kinds of correlation functions of 
the pair annihilation and coalescence models are related to one another in any dimension. 



This method also reproduces all results in Ref. [22| regarding the equivalence of stochastic 
systems by adjusting fi and also finds the initial condition relation. 

Next, we apply the generating function method to find the equivalence between three- 
particle annihilation models. To our knowledge, this equivalence has not been studied in 
the literature, although the full renormalization-group study of this model can be found in 
Ref. ^. The dynamics of the three-particle annihilation model is summarized in Table ITTl 
The procedure to find the equivalence relations is the same as those explained above. First, 
we find the Hamiltonian and the corresponding evolution operator. The corresponding 
Hamiltoman is = Eti [h^, + K', + + ^^.) with 

= D [pM - Pri+eJ - OrioUe, + (1 " Pn)pn+e, " aj,a„+ej , (33a) 
^n,i ~ ['^PnPn+eiPn+2ei ~ 0,nPn+eiPn+2ei ~ PnPn+eiO'n+2ei ~ PnO'n+eiPn+2e.i\ ; (33b) 

^ri,i ~ ~^ \^PnPn+eiPn+2ei ~ 0'nO'n+e.iPn+2ei ~ (^nPn+eiO'n+2ei '~ Pn(^n+ei(^n+2ei\ ; (33c) 
^n,i ~ ^3[PnPn+eiPn+2ei — 0'nO'n+eiO'n+2e J, (33d) 

and the evolution operator is £ = Ylm Ylii ^n,i with 

f d d \ d d 



1 _ _ _ ^ d d d 



1 ____ _ _ d d d 

+ o (2Al + 3A2 + 'i\'i){^nVn^e.i + V^n+e, V^ri+2e, + ^n+2e.i^nj 



(34) 



3 dipn Oifn+e, OVn+2e, 

, , . _ _ d d d 

+ (,Ai + A2 + A3jV5riV5ri+e,V5ri+2e,77— TTT TTT • 

By rescaling ip = ^ip', we find the parameter relations 

D = D,~Xi + 2X2 + 3A3 = p'\Xi + 2X2 + 3A3), 

2Ai + 3A2 + 3A3 = p-\2Xi + 3A2 + 3A3), (35) 

Ai + A2 + A3 = Ai + A2 + A3. 

First consider 3^4 only (Ai = A2 = 0). Unfortunately, the only possible value of yU is 1, 
which means there is no equivalent procedure to 3A 0. However, if we set A2 = A3 = 0, 
we find 

Ai = Ai (3 - ^) , A2 = Ai^i±i[^^, A3 = Ai(l - /i-l)^ (36) 
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where the vahd range of/iis2/3<;U<l. 

Up to now, we have found the equivalence for single-species react ion- diffusion systems, 
but the formalism explained above can be extended to multispecies problems. As an exam- 
ple, let us consider a rf- dimensional two-species annihilation model. There are two kinds of 
species, A and B. Both particles have the same diffusion constant D. When two different 
species form a pair, one of the following reactions may occur. Both particles are removed, 
or only one particle is annihilated. For the-one particle annihilation process, there are two 
possibilities. The remaining species stays where it was or moves to the other site. For these 
dynamics, we assign transition rates A (pair annihilation), rj (one-particle annihilation with- 
out location change), and ( (one-particle annihilation with location change). The dynamics 
is summarized in Table UTTl 

The Hamiltonian for this model is H = X]f=i Hn,i with 

Hn,i ~ ^{A-nBn+ei '~ O-nPn+ei + B^An+ei ~ bnO-n+ei): (37b) 
Hn,i ~ Vi^nBn+Si ~ 0,nBn+ei + A^Bn+et ~ A^bn+ei (37c) 
-\-ByiAyij^g^^ bnAn-^e- -\- ByiA^j^^^ -Sj^fl^_(_e^ ) , 

~ C{AnBn+ei ~ b^a^Jjn+ei + A^Bn+e^ — CLnO!^ <^9n+eibn+e.i (37d) 
tt^bnOin+ei ~l" B^An+ei brib^j^g^^n+ei) } 

where An {Bn) is the number operator for species A (B), and Vn is the vacuum projection 
operator with the form {)„ = 1 — A^ — Bn- We use a different symbol for the different 
species and drop the species indices. The commutation relations for the operators in Eq. 
()37|) correspond to the M = 2 case of Eq. (jlll) . The evolution operator takes the form 
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^= Er.E^tl'^r.,^ with 



d 



d 



d d 



d d 



+ (D-r^) 

+ (C + 2ri + X-dD) 



d d 



d d 



diPn dl/^n+ei dl/j^ dip. 



n+ei 



— d d — d d 



di/jn dip. 



n+ei 



n+ei 



— 00— 00 

{<fn + i^n+e,) ^—^7 + {^n + (Pn+e,) 



Oipn Oi/jn+e 



Olpn dipn+ei 

( _ - - _ 

+ (2C + 277 + A - 2(iL>) </^„V'n+ei;^— ^7 V i^n^n+e,^^-^ 

V Oipn Oi>n+ei O^J^i OiPn+. 

where (^(■0) corresponds to species A (B). By rescaling, we find the relations 

D^D, D-fj^-{D-rj), 
A* 

C + 2fi+~X-D=-{C + 27] + X-D), 
A* 

2C + 2?) + A - 2L' = 2C + 277 + A - 2D, 



or 



A* 

A = A+ -{C + ri + X), ri = - + - D. 



If we set ?7 = C = 0, we can find 



(38) 



(39) 



(40) 



(41) 



li li H 

If A > D, we obtain equivalent stochastic processes with the pair annihilation model A-\-B — > 
for 1 < < 2. Other equivalent relations can be found by adjusting ^ for different values 
of untilded parameters. 



IV. SUMMARY AND DISCUSSION 



In summary, we have presented the path integral representation of the generating func- 
tion for classical stochastic exclusive systems. Using this formalism, we have shown the 
equivalences of some stochastic systems of single- and two-species models. Furthermore, 
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since our formalism is insensitive to the choice of the function / in Eq. ()16|) . we can have 
different action forms for one system depending on the choice of /. Although they may 
have different forms, the long-time behavior of physical quantities derived from each action 
form should be the same. For example, if we use the prescription Eq. (jl8c|) in Eq. 
the action is equivalent to the action with the choice of /(<^) = <^ in Eq. and with the 
change of fields, such as <^ ^— e*^ and (p ^— e""^!^, performed. In fact, if we use the function 
/((^) = e'^ and expand the function / in terms of (p keeping only terms up to ip^ order, 



we obtain the same action as the one derived in Ref. [17[. In many cases, any term which 
contains <^ higher than second order is irrelevant. The method introduced in Ref. is 
legitimized by the force of the renormalization-group argument. In some cases, usually in 
react ion- diffusion systems, the canonical dimension of is 0, 3|, so it is not clear whether 
we can neglect the higher-order terms after expanding /. In this situation, one should use 
Eq. p8a|) or Eq. ()18b|) to ensure the completeness of the formalism. As an example, let us 
consider again the model described by the reaction given in TablelHwith r/ = 5 = a = 0. If 
we use Eq. fllSbjl . the action becomes (up to relevant terms) 

S = j dt(fx y{dt- DV^) up + 2A^(/?2 + \fp?} , (42) 



where we neglect the term D 



■ El 



'yp>f' due to its irrelevance. Equation (jl^ is exactly the 



same as the action in Ref. obtained using a bosonic formalism. Since we expect that 
boson and hard-core boson annihilation models should belong to the same universality class, 
this coincidence is not surprising. If we perform the change of fields explained above, and 
keeping terms up to 0[}p^\ the action takes the form 

S' = \ dtd'^x{<p{dt - DV^)^ - DipiVipf 

J (43) 

+ 2Xipip^ - 2A(^V^}, 

which can also be found by the method introduced in Ref. 17[. These two actions, Eq. 

and Eq. PHj). look different, although they describe the same system. This is due to 
the dropping of higher-order terms of ip. To describe the system properly, we have to keep 
all the higher-order terms of ip which are marginal. Then, two actions will give the same 
long-time behavior of physical quantities. 

Although we have concentrated on react ion- diffusion systems, the applicability of this 
forma 
DLG 



ism is not restricted to reaction-diffusion systems. We found Langevin equations for 
and other lattice-gas- type models. This will be published elsewhere j^- 



13 



Acknowledgments 

This research was supported by the Cathohc University of Korea research fund 2004 and 
Grant No. R05-2003-000- 12072-0 from the BRP program of the KOSEF. 



[1] B. Chopard and M. Droz, Cellular Automata Modeling of Physical Systems (Cambridge Uni- 
versity Press, Cambridge, 1998). 

[2] J. Marro and R. Dickman, Nonequilibrium Phase Transitions in Lattice Models (Cambridge 
University Press, Cambridge, 1999). 

[3] M. Doi, J. Phys. A 9, 1465 (1976). 

[4] P. Grassberger and M. Scheunert, Fortschr. Phys. 28, 547 (1980). 

[5] L. Peliti, J. Phys. (Paris) 46, 1469 (1985). 

[6] B.P. Lee, J. Phys. A 27, 2633 (1994). 

[7] B.P. Lee and J. Cardy J. Stat. Phys. 80, 971 (1995). 

[8] M.W. Deem and J.-M. Park, Phys. Rev. E 57, 2681 (1998); J.-M. Park and M.W. Deem, ibid. 
57 3618 (1998). 

[9] R. Dickman and R. Vidigal, J. Phys. A 35, 7269 (2002); Braz. J. Phys. 33, 73 (2003). 
[10] S.A. Janowsky Phys. Rev. E 51, 1858 (1995). 

[11] I. Ispolatov, P.L. Krapivsky and S. Redner, Phys. Rev. E 52, 2540 (1995). 
[12] K.E. Bassler and D.A. Browne, Phys. Rev. E 55, 5225 (1997). 

[13] R.K.P. Zia and B. Schmittmann, Phase Transitions and Critical Phenomena, edited by C. 

Domb and J.L. Lebowitz (Academic Press, New York, 1995), Vol. 17. 
[14] V. Brunei, K. Oerding, and F. van Wijland, J. Phys. A 33, 1085 (2000). 
[15] P.-A. Bares and M. Mobiha, Phys. Rev. E 59, 1996 (1999). 
[16] F. van Wijland, Phys. Rev. E 63, 022101 (2001). 
[17] S.-C. Park, D. Kim, and J.-M. Park, Phys. Rev. E 62, 7642 (2000). 

[18] S.-C. Park, D. Kim, and J.-M. Park, Phys. Rev. E 65, 015102(R) (2002); S.-C. Park, J.-M. 

Park, and D. Kim, ibid. 65, 036108 (2002); S.-C. Park and J.-M. Park, ibid. 67, 010103 (2003). 
[19] C.W. Gardiner, Handbook of Stochastic Methods, 2nd cd. ( Springer- Verlag, Berlin, 1990). 
[20] H. Risken, The Fokker-Planck Equation (Springer- Verlag, Berlin, 1984). 



14 



[21] N.G. van Kampen, Stochastic Processes in Physics and Chemistry, enlarged ed. (Elsevier, 
Amsterdam, 1997). 

[22] M. Henkel, E. Orlandini, and G.M. Schiitz, J. Phys. A 28, 6335 (1995). 
[23] T. E. Harris, Ann. Prob. 2, 969 (1974). 
[24] H. Hinrichsen, Adv. Phys. 49, 815 (2000). 
[25] S.-C. Park and J.-M. Park (unpublished). 

TABLE I: Reaction-diffusion processes of single species and their rates. 

Diffusion ^ $A with rate D 

Pair annihilation AA 00 with rate A 

Coalescence A A A$ with rate rj/2 

Coalescence A A with rate rj/2 

Death ^ ^ with rate S 

Branching $A — >■ A A with rate a/2 

Branching ^0 AA with rate a/2 



TABLE II: Reaction-diffusion processes of the three-particle annihilation model and their rates. 

A$ ^ (l)A with rate D 

AAA AM, MA, or AA$ with rate Ai/3 

AAA ^00, or 00A with rate A2/3 

AAA 000 with rate A3 
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TABLE III: Dynamics of the two-species annihilation model and respective rates. 



4(/l 


Willi Idle J-/ 


R(7I (TIR 


Willi Idle J-/ 


4 R _v p)(7i 


Willi Idle /\ 


BA 00 


with rate A 


AB yl0 or 05 


with rate rj 


BA 50 or 0^ 


with rate rj 


^45 ^ 50 or M 


with rate ^ 


BA A0 or 05 


with rate ^ 
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